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Abstract 



We consider Cauchy problem for Fourier transformation of 3-dimensional Navier-Stokes system with 
■ zero external force. Using initial data purposed by Dong Li and Ya. G.Sinai in [I] we implement self- 

similar regime producing fast growing behavior of the energy of solution while time tends to critical 
value T rr . 



< 

1 Theoretical results used. 

Fourier transform of the 3D Navier-Stokes system with purely imaginary initial data can be written 
in the form 

00 ! v(k, t) = e- tW2 v(k, 0) + / e -(*- s )l fc l 2 / (k, v(k - I, s))P k (v(l, s))dlds (1) 

where k £ M 3 corresponds to Fourier mode of the solution, t £ [0, T cr ] C K+ is time, v(k, t) corresponds 



o 



to imaginary part of the solution, (•, •) denotes Euclidian inner product. P k is a Leray projection on 



the subspace orthogonal to k and has the form P k = Id rrrr^- Viscosity supposed to be equal 1, 



external forcing is put to zero. Incompressibility condition means that for all t G [0, T cr ] 



^! v(k,t)±k forany/teM 3 (2) 

For given function co(k) consider one-parametric family of initial conditions VA,(k,0) = Aco(k). 
Accordingly to Ya.G. Sinai [2], [3] for any cq with fast decay at infinity, for every A there exists 
such constant T cr that for t £ [0, T cr ] there exists unique solution of the system ([T]) which can be 
represented as a power series of the parameter A: 



va 



(k, t) = Ae- m \ {k) + f e -(*- s )l fc l 2 A p h p (k, s)ds (3) 

n P>1 



Substituting ([3]) into ([T]) one can express coefficient h p (k,s) trough coefficients with fewer indices 
and from initial data cq with the following recurrent relations: 
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h 2 (k, s) = J (I, c (k - l))P k c (l) exp{-s{\l\ z + \k - l\ z }dl 

s 

h p (k,s)= J dsi J (l,co(k — l))Pkh p _i(l, si) exp{—(s — si)\l\ 2 — s\k — l\ 2 }dl+ 

R3 

s 

+ J ds 2 j \l,h p -i(k - l,s 2 ))PkCo(l)exp{-s\l\ 2 - (s - s 2 )\k - l\ 2 }dl+ 

R3 

s s 

+ Yl / dsi / ds2 / {^ h Pi{ k -hsi))P k hp 2 {l,s 2 ) (4) 

Pi+«0 R3 

exp{— \k — l\ 2 (s — s\) — \l\ 2 (s — s 2 )}dl 

Consider initial data with compact support C located in the ball of radius r with the center at 
the point (0, 0, R). From recurrent relations (4) it follows that the coefficient h p is supported in 
the algebraic sum C + C + ■ ■ ■ + C. Then one can fix 1 < r such that for the ratio large enough 

p times 

supports of the coefficients h p and h„/ p' ^ p does not intersect. In more details, support of each 
function h p lies in the ball of radius 0(l)^/p with center at the point (0, 0, R ■ p). 

We reproduce sketch of the arguments of Ya.G. Sinai and Dong Li [lj. Third component of 
coefficient h p can be expressed through first and second ones due to the incompressibility condition 
([2|). Since all terms in (4) where either p\ or p 2 is small comparable to ^/p can be estimated with 
some function which decays with p faster than exponent one can concider only the part where both 
Pi ~ VP an d p 2 ~ y/p- Change variables according to the ball containing the support of function h p 
and consider a function g p (y, s) = h p ((0, 0, R) + VRy, s). After renormalisation, we get 

, . . , x )\ — 

P1+P2=P v\v 2 j v V2 \pi ■ y pi 

P1,V2>VP R (5) 

+{V2 - x 2 )g$((y - x). — ,s) + x 1 g p 1) ({y - x). —,s) + x 2 gf^{(y - x). — ,s) ) dx 

V pi V pi y pi J 

Write down functions g p (y,s) in the form 

g p (y, s) = A p (s)pa 1 exp \-^-(\yi\ 2 + M 2 )} sfoi exp {o- 2 |y 3 | 2 } 

x (6) 

Hi{y) + S^iy, s), H 2 (y) + #\y, s), #\y, s) 



9p(y^ s ) = ^ 2 V] -f"2 [ 9 P2 (, x x P^,s) (— (yi - x\)g$((y 
R „,irf_„ P1P2 J V P2 \P1 



While p tends to infinity one can show that functions 5^f\y,s), 8^ \y,s) and 6^\y,s) became 

Pi 

o(l). Sum in the right-hand side of the ([5]) can be considered as integral sum with variable 7 = — . 

P 
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For the components Hi and H 2 we get from dSJ),© (for details see [U Eq.21])) 



2vr |c7 (1) |x| 2 
(H 1 (x),H 2 (x)) = — exp 




(1)^2 



R2 



(2vr) 2 7 (l-7) 



exp 



a (l)| X _y|2 a (l)| y |2 



2 7 



2(1-7) 



-(1 - 7)» ^^i(^) + ^^tf 2 (^) ) + 
1 V7 1 



+(7(1 - ( intfi(^) + y 2 ^ 2 (^) 



(7) 



Ya.G. Sinai and Dong Li in [lj show that operator acting on functions {Hi,H 2 } given by right-hand 
side of the equation (|7|) has 6 unstable eigen values and 4 neutral. In other words for any T cr > 
there exist such constant A(T cr ) and 10-parametric family of initial data co(k) such that corresponding 
functions g p after renormalisation tends to some limit form g(y,s), while p tends to infinity. Energy 
of the solution of the equation (JTJ), constructed from functions g(y, s) tends to infinity as t —* T cr . In 
PQ there was obtained an asymptotics for energy growth 



8(t) 



\v(k,t)\ 2 dk 



(T cr - ty 



Proof of the existence theorem in [T] does not contain particular form of the initial data and corre- 
sponding solution. Thus numerical study of such solutions is of known interest. 



2 Computational scheme and main result. 

Differentiation by time t in the lefthand side of the equation (pQ) represents traditional form of the 
ordinary differential equation for the vector-valued function v(k,t) for any fixed k. Convolution-type 
operator in the righthand side has non-local nature i.e. value of solution v(k,t) in particular point 
(k,t) depends on values of the solution in each point k G M 3 . Computational model is constructed 
in such domain f2 C M 3 that for all ^ t < T cr values of function v(k, t) outside of the Q can be 
neglected. For the discretisation of the arguments t and k we consider difference operator for t and 
computation with discrete argument instead of integration in the righthandside of ([T|) 
Rewrite space integration term in the equation (pQ) using incompressibility condition 

J (k, v(k - I, s))P k v(l, s)dl = P k j (I, v(k - I, s))v(l, s)dl 

R 3 R 3 

Expression in the righthandside can be treated in terms of convolutions of three-dimensional arrays. 

Look for solutions of the equation ([1]) with initial data located in the ball of radius 0(y/R) centered 
in the point ko = (0, 0, R) and having in this support form of linear combination of rescaled Hermit 
functions. In more details we compute initial values v(k,0) = (v 1 (k) , v 2 (k) , (k)) in points of three- 
dimensional grid with given accuracy using expression 

^)=nEv eW (^#). ^ =i ' 2 > 3 ( § ) 

i=l m=l \ V / 



Here k l , v J correspond to components of the vectors k and v respectively, He^ m > f ^ j is a 



Hermite function of order m with argument 3 ^ k ^ ^ and {Xij m } corresponds to the set of coefficients. 
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Available computational resources made it possible to consider grids of order (const • 10 2 ) x (const • 
10 2 ) x (const • 10 3 ), where constants may vary from 1 to 2. According to this accuracy of grid in each 
coordinate vary from 0.1 to 0.2 respectively. For the case R == 5 support of the solution v(k, t), which 
we were able to observe lie in the domain [—15, 15] x [—15, 15] x [1, 200]. 

For calculation of three-dimensional convolutions we used Intel freeware library for non-commercial 
using Math Kernel Library. 

Since due to [2] for the solution v(k,t) we have asymptotic one can estimate accuracy of 

the computational scheme. With space resolution near 0.2 relative accuracy can be estimated with 
10%. 

Tornado— type solutions. 

Main result of this paper consists in presenting particular form of the initial conditions 
for which energy of the solution £{t) has poly normally fast growthrate while t tends to 

j-cr- 

Rough, due to our computational resources, estimate of the order of growth rate of the energy 
of the solution is far from the teoretical l/(T cr — i) 5 . Observed asymptotics for this order with T cr 
defined with respect to accuracy of our calculations are near 15 - 20. Thus one may suppose that 
formation of blowing up solution, consturted from the functions g(y, s) as it proposed in pQ is not 
unique. 

For the moment we have no examples of the solutions with growth rate ~ 5. For large set of 
variations of parameters D, Aij m and form of initial condition solution decay in time. 

However, founded behavior of numerical solution is stable for slightly perturbed (near 25%) dis- 
cretisation step in time, space and form of initial data. Form of the blowing solution for such perturbed 
initial data differs very slightly, and time of the solution transforms with natural rescaling of (pQ). 



3 Results. 



Below we present as illustration computed results for the support located in the ball of radius R == 5, 
centered at the point (0,0,5). Maximal order of Hermite polynom is taken to be D == 3, starting 
discretisation time step is dt == 0.001, space resolution is 0.14907 thus number of points of discretisation 
in X and Y coordinates are nX == nY == 160, number of discretisation points on third coordinate is 
nZ = 1440. We present square root of the integrated square of the absolute value of the solution 
(mynorm == y/£) in consequent time moments. From t == 0.000 to t == 0.042 energy monotonically 
decreased. 

t=0.000 A4(i)=1.00000 t=0.042 A*(t)=0.3664 t=0.045 .M(i)=19.972 
t=0.040 .M(i)=0.32191 t=0.043 A*(t)=0.9500 t=0.046 A4(i)=101.20 
t=0.041 7W(t)=0.31746 t=0.044 7W(t)=4.1879 t=0.047 M(t)=698.57 
In the second table we present results of the computations with the same initial data with different 
discretisation parameters. Space discretisation here is 0.20328, time step dt == 0.0008, number of points 
in X and Y coordinartes are nX == nY == 120 and for third coordinate nZ = 1080. 



t=0.0400 
t=0.0408 
t=0.0416 



M(t)- 

M(ty- 
M(ty 



:0.3226 
=0.3165 

:0.3105 



=0.0424 
=0.0432 
=0.0440 
=0.0448 



Mity 
M{ty 
Mity 
M(ty 



=0.3047 
=0.2996 
=0.3017 
=0.3982 



=0.0456 
=0.0464 
=0.0472 
=0.0480 



M(t)-- 
M(t)-- 

M(ty 
M(ty 



=1.1492 
=4.6551 
=20.669 
=1027.7 
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Figure 1: Absolute Value in uniform and log scales 
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Figure 4: First cloud. Main contribution. 




Figure 5: Second and Third clouds. 
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